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Abstract 

Two-particle correlations have been measured for identified ir~ from central 158 AGeV 
Pb+Pb collisions and fitted radii of about 7 fm in all dimensions have been obtained. A multi- 
dimensional study of the radii as a function of kr is presented, including a full correction for 
the resolution effects of the apparatus. The cross term Ro Ut _i ong of the standard fit in the 
Longitudinally CoMoving System (LCMS) and the vl parameter of the generalised Yano- 
Koonin fit are compatible with 0, suggesting that the source undergoes a boost invariant 
expansion. 

The shapes of the correlation functions in Qi nv and Q spa ce = y ^Q x + Qy + Q 2 z have been 
analyzed in detail. They are not Gaussian but better represented by exponentials. As a 
consequence, fitting Gaussians to these correlation functions may produce different radii 
depending on the acceptance of the experimental setup used for the measurement. 



1 Introduction 

The study of Bose-Einstein correlations between pairs of identical hadrons is an essential tool 
to obtain information on the space-time evolution of the extended hadron sources created 
in heavy ion collisions In particular, a strong correlation between the momenta and 
the space-time production points of the particles suggests expanding sources as predicted 
by hydrodynamic models . The dynamical evolution of such systems can then be studied 
with interferometry via selection on the transverse momenta and rapidity of the correlated 
particle pairs. 

In this paper, we present the analysis of two-particle correlations of identified tt~ measured 
in the WA98 experiment for central 158 AGeV Pb+Pb collisions at the CERN SPS. 



2 Experimental setup and data processing 

The WA98 experiment shown in Fig. |l| combined large acceptance photon detectors with a 
two arm charged particle tracking spectrometer. The incident 158 AGeV Pb beam interacted 
with a Pb target near the entrance of a large dipole magnet. Non-interacting beam nuclei, 
or beam fragments were detected in a forward calorimeter located at zero degree. A mid- 
rapidity calorimeter measured the total transverse energy in the rapidity region 3.2 < rj 
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< 5.4, which was also used in the trigger for online centrality selection. The Plastic-Ball 
calorimeter measured the fragmentation of the target, and silicon detectors were used to 
measure the charged particle multiplicity. The photon detectors consisted of a large area 
photon multiplicity detector and a high granularity lead-glass calorimeter for single photon, 
7T°, and 7] physics Q. 

The charged particle spectrometer made use of a 1.6 Tin dipole magnet with a 2.4x 1.6 m 2 
air gap for magnetic deflection of the charged particles in the horizontal plane. The results 
presented in this paper are taken from the 1995 WA98 data set obtained with the negative 
particle tracking arm of the charged particle spectrometer. The second tracking arm was 
added to the spectrometer in 1996 to measure positive particles 0. The first tracking arm 
consisted of six multistep avalanche chambers with optical readout [f| located downstream of 
the magnet. The active area of the first chamber was 1.2x0.8 m 2 , while that of the other five 
was 1.6x1.2 m 2 . The chambers contained a photoemissive vapour (TEA) which produced 
UV photons along the path of traversal of the charged particles. These were converted into 
visible light via wavelength shifter plates. On exit the light was reflected by mirrors at 45° to 
CCD cameras equipped with two image intensifiers. Each pixel of a CCD viewed a 3.1x3.1 
mm 2 area of a chamber. In addition, a 4x1.9 m 2 Time of Flight wall positioned behind the 
chambers at a distance of 16.5 m from the target allowed for particle identification with a 
time resolution better than 120 ps. 

Fig. ^ shows the Monte Carlo generated p^-rapidity acceptance for 7r~. The acceptance 
ranges from y=2.1 to 3.1 with an average at 2.70. The momentum resolution of the spec- 
trometer was Ap/p=0.005 at p=1.5 GeV/c, resulting in an average accuracy better than or 
equal to 10 MeV/c for all the Q variables used in the interferometry analysis and defined 
in section 5: a(Q inv )=7 MeV/c, a(Q T o)=W MeV/c, a(Q T s)=5 MeV/c, a(Q L )=3 MeV/c, 
a(Q T )=8 MeV/c, a(Q )=5 MeV/c. 

The analysis of the complete 1995 data set is presented here. These data have been taken 
with the most central triggers corresponding to about 10% of the minimum bias cross section 
of 6190 mb. Severe track quality cuts were applied at the expense of statistics resulting in 
final samples of 4.2 xlO 6 7r~ for the correlation analysis and 4.6xl0 5 tt~ for the single particle 
spectrum. 

3 Single particle spectra 

The mr=y / m 2 + p\ distribution of identified tt~ , averaged over the rapidity acceptance, is 
shown in Fig. |3[ The data were corrected for geometrical acceptance and efficiency of the 
chamber-camera-Time of Flight system using a full simulation of the experimental setup. The 
parameters of the simulation were optimized in an iterative way by comparing various distri- 
butions with the real data. The simulated data were then treated exactly like the real data. 
The measured 1 / mj- dN j dm? dy distribution was then fitted to the form Cexp(— m^/T), 
expected for a source in thermal equilibrium ||. Such fits were applied to the data for dif- 
ferent ranges of rriT, such as the one shown in Fig. [3|. These fits do not reproduce the overall 
concave shape of the data, which is partly due to particles originating in resonance decays 
and could also be an indication of transverse flow (7j. The shape of the ir~ tut distribution 
was found to be in good agreement with that of 7r° obtained in the lead-glass calorimeter || . 
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4 One-dimensional interferometry analysis 



For the Bose-Einstein correlation studies, the data were Coulomb corrected in an iterative way 
|8). The Gamow correction was abandoned as it overcorrects the data for Qi nv in the range 
of 0.1 to 0.3 GeV/c. A fit of the form 1 + Aexp[-Q? w i?? n J was made to the Q% nv correlation 
function yielding Rinv = 6.83 ± 0.10 fm and A = 0.307 ± 0.008. An expanded view of the 
correlation distribution (Fig. ||) shows that the Gaussian fit used (full line) is not perfect, 
especially in the Qi nv range of 40 to 80 MeV/c where the tail of the experimental distribution 
shows an excess which is not well reproduced by the fit. In addition to this Gaussian fit made 
over the whole range of Qi nv , Fig. ^| shows also different Gaussian fits using data in the ranges 
25 MeV/c < Q inv < 200 MeV/c (dashed line) and 40 MeV/c < Q inv < 200 MeV/c (dotted 
line). These fits do not coincide. Different radii are then obtained for different starting points 
of the fit because the shape of the distribution is not Gaussian. This effect is independent 
of the severity of the track selection, and is therefore not due to spurious tracks. This is 
summarized in Fig. |6] where Ri nv and the corresponding A are plotted as a function of the 
lower bound of the fit. There is a statistically significant drop when using a Gaussian fit. 
A similar behaviour is observed when, instead of Qi nv , Q space = \J Q"i + Qy + Ql is used, 
calculated in the longitudinally comoving system (LCMS) and fitted with 1+Aexp[-Q|i?|]. 
This method of fitting in varying ranges has a good sensitivity to the shape. It has been 
repeated by replacing the Gaussian fit by an exponential fit of the form l+X e eyip[-2Qi nv R e ] 
where the factor 2 is added to make the radius R e more comparable with Ri nv . The results 
(Fig. 0) show that the stability is better with the exponential fit. Fig. || directly compares the 
Gaussian and exponential fits for Qi nv . Although the Gaussian fit still gives an acceptable 
% 2 /d.o.f., the exponential fit is better everywhere. A similar conclusion is reached when the 
first data point is excluded from the fit. This result is not based on the first bins which might 
be more affected by systematics due to large Coulomb correction or noise correlated with 
the true track signals in the chambers. It is rather based on the high statistics tail of the 
distribution which contributes in a different way to a Gaussian or an exponential fit. This 
quasi exponential behaviour is expected by different models including resonance decays ||. 
As a consequence small acceptance experiments may obtain a larger radius if a Gaussian fit is 
used because they are less sensitive to the tail. On the contrary large acceptance experiments 
have higher statistics at large Q- values, and the Gaussian fit will yield lower values of the 
radius. 

5 Multi-dimensional interferometry analysis 

The multi-dimensional analysis has been done with Gaussian fits to allow comparison with 
other experiments. Two different parameterizations have been used in the LCMS: 

a) The standard fit in the 3-dimensional space of momentum differences Qts (perpen- 
dicular to the beam axis and to the transverse momentum of the pair), Qto (perpendicular 
to the beam axis and parallel to the transverse momentum of the pair), and Ql (parallel to 
the beam axis) |pXf |. The fitted formula 

C 2 = 1 + Aexp[-Qy 5 i?y 5 - Qto r to ~ Ql r l ~ ^QToQlR-lut-long] 
includes a cross term in QtoQl as predicted (IHl . 
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b) The generalized Yano-Koonin (GYK) fit 12] in the Qo (energy difference of the pair), 
Qt-,Ql space according to 



C 2 = 1 + XbxpI-QtRt + (Qo - Ql)Rl - (Q-U) 2 (R 2 o + R\)] 

where U = 7(1, 0, 0, vl), 7 = 1 / \J 1 — v\ with vl in units of c=l. 

In the GYK approach, the radius parameters remain invariant under longitudinal Lorentz 
boost, the parameter vl connecting the "arbitrary" measurement frame (LCMS) to the Yano- 
Koonin frame. In addition, the extraction of the duration of emission, Rq, is straightforward. 

The consequence of the finite resolution in the measurement of the Q variables is an 
underestimate of the radii and A parameters. Morever, as the resolution is different for each 
Q variable, this causes a bias which varies from parameter to parameter, leading to errors in 
the interpretation of the results in a multi-dimensional analysis. It is therefore essential to 
take into account the effect of the resolution in the fitting procedure. One way to do this is 
to replace the formula C2{Q) used to fit the data by 

C r 2 c (Q) =JJJ r(Q,Q>) C 2 (Q') dQ> 

which is the convolution of C2(Q) with the resolution function r(Q,Q'). The resolution 
function is chosen to be Gaussian: 

r(Q, Q>) = l/(2vr) 3 / 2 1/| VI 1 / 2 exp[-l/2 (Q - Q'f V~ l (Q - Q')} 

The diagonal elements of the covariance matrix V are equal to the square of the resolution of 
the different Q variables and are estimated separately as a function of = |pri+PT2|/2 of the 
pairs. The non-diagonal elements are neglected. For the one-dimensional Gaussian fit case 
with Q = Qi nv , the resolution corrected values of the fitted parameters are Ri nv = 7.30 ±0.12 
fm and A = 0.328 ± 0.009. 



Table 1: 3- dimensional analysis 



Standard fit in LCMS 


Generalized Yano-Koonin fit 


Rts 


= 6.41±0.13 fm 


Rx 


= 6.54±0.11 fm 


Rto 


= 6.60±0.16 fm 


R 


= 0.01±0.69 fm 


Rl 


= 7.50±0.18 fm 


R4 


= 7.51±0.18 fm 


A 


= 0.350±0.010 


A 


= 0.325±0.009 


-^out—long 


-1.0±1.3 fm 2 


vl 


= 0.03±0.05 


xVd.o.f. 


= 1.06 


X 2 /d.o.f. 


= 1.02 



The results of the multi-dimensional fits are presented in Table |T| for the full 1995 data 
sample. A multi-dimensional analysis as a function of /ct, both with the standard 5-parameter 
fit and with the GYK fit is shown in Figs. ||, |9[ and |l(]. 

The Rts and Rl parameters from the standard fit are found to be compatible respectively 
with Rt and R4 from the GYK fit. The cross term R 2 out i ona from the standard fit and vl 
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from the GYK fit are compatible with 0. In a source undergoing a boost invariant expansion, 
the local rest frame coincides with the LCMS. Both the cross term and vl estimated in the 
LCMS are then expected to vanish [12|. As this is the case, it suggests that the source 
seen within the acceptance of the experiment has a strictly boost invariant expansion. The 
strong decrease of the longitudinal radius Rl or R± with kr compared to the behaviour 
of the transverse radii Rt, Rts, Rto suggests a longitudinal expansion larger than the 
lateral expansion. The longitudinal radius Rl is shown with a fit of the form 1/^/mj with 



mT=ym,2 + fcj. inspired by the hydrodynamical expansion model. Using Rl = Toy/To/rriT 
with a freeze out temperature To of 120-170 MeV/c, we may extract a freeze out time tq in 
the range of 7.5-8.9 fm/c. Finally, the Rq parameter from the GYK fit, which reflects the 
duration of emission, is compatible with for all kx bins, excluding a long-lived intermediate 
phase. 

Two other experiments, NA49 and NA44, have studied charged particle interferometry 
in Pb+Pb collisions at CERN energies. The WA98 analysis is in good agreement with the 
NA49 results [13 1, when the comparison is made for the same mean y range of 2.70, although 
WA98 has used identified ir~ while the NA49 analysis used unidentified negative particles. 
Only the Ro parameter tends to be smaller in WA98. The direct comparison with the NA44 
experiment is not possible because NA44 and WA98 do not have the same y range. The 
smaller radii measured by NA44[14] can be explained by the larger y range of its acceptance 
(3.1<y<4.1). 



6 Conclusion 

The analysis of the two-particle correlation of identified 7r~ from central Pb+Pb collisions at 
158 AGeV gives fitted radii of about 7 fm. This should be compared to the equivalent rms 
radius of the initial Pb nucleus of 3.2 fm, which indicates a large final state emission volume. 
The one-dimensional correlation functions analyzed in terms of Qmv or Q space 

are not 

Gaussian. They are better represented by exponentials. This study is based on the tail of 
the distributions and not on the first bins which might be subject to systematic effects. One 
possible explanation is that this behaviour is due to resonance effects. Fitting Gaussians to 
these correlation functions may produce different results depending on the acceptance of the 
experimental setup. 

The generalized Yano-Koonin analysis gives similar results to within the error bars as the 
standard 3-dimensional analysis in the LCMS. 

The cross term R^ mt _i ong is found to be compatible with in the LCMS and the same 
is true of vl in the GYK fit. This suggests that the source undergoes a boost invariant 
expansion. 

A clear dependence of the longitudinal radius parameter on kx is observed, suggesting a 
larger longitudinal than transverse expansion of the source. In addition the short duration 
of emission disfavours any long-lived intermediate phase. 
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WA98 Experimental Setup 
158 A GeVPb+Pb Collisions 
at the CERN SPS 
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Figure 1: Experimental setup. 
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Figure 2: p T -rapidity acceptance. 
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Figure 3: ttit distribution (the fit is explained in the text). 
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0.02 0.04 0.06 0.08 0.1 



Q inv [GeV/c] 

Figure 4: Q inv distribution with Gaussian fits. The full line is a fit made with all the points. The 
dashed line is a fit for Q inv > 25 MeV/c and the dotted line is a fit for Q inv > 40 MeV/c. 
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\ X 2 /dof=0.89 


X 2 /dof = 2.12 
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QinvlGeV/c] 

Figure 5: Q inv distribution with the Gaussian fit (full line) and the exponential fit (dashed line). 
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Figure 6: Gaussian fit stability. 
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Figure 7: Exponential fit stability. 
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Figure 8: 5-parameter fit as a function of kr- The horizontal bars indicate the bin width and the 
points are plotted at the average kr of the bin. The dashed line in the bottom left figure shows 
the result of a fit explained in the text. 
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Figure 9: Comparison of the Rts and Rl parameters from the standard fit respectively with R T 
and i? 4 from the GYK fit as a function of k T . 
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Figure 10: The cross term R 2 out _i ong from the standard fit and Vl from the GYK fit, both estimated 
in the LCMS, are compatible with for all k T bins. 
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